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Abstract 

We study a self-interacting scalar field theory in the presence of a (5-function background poten- 
tial. The role of surface interactions in obtaining a renormalizable theory is stressed and demon- 
strated by a two-loop calculation. The necessary counterterms are evaluated by adopting dimen- 
sional regularization and the background field method. We also calculate the effective potential 
for a complex scalar field in a non-simply connected spacetime in the presence of a (5-function 
potential. The effective potential is evaluated as a function of an arbitrary phase factor associated 
with the choice of boundary conditions in the non-simply connected spacetime. We obtain asymp- 
totic expansions of the results for both large and small 5-function strengths, and stress how the 
non-analytic nature of the small strength result vitiates any analysis based on standard weak field 
perturbation theory. 
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I. INTRODUCTION 



A. Dedication 

It is a great pleasure to contribute this article to honour Stuart Dowker's many achieve- 
ments in theoretical physics. I first knew Stuart by reputation long before I ever had the 
chance to meet him. When I was a postgraduate student my main interest was in quantum 
field theory in curved spacetime. Long before the internet and the arXiv preprints in high 
and theoretical energy physics used to be sent by post to SLAC where they would appear 
on a list that was sent out weekly to subscribers. At that time there were four main groups 
in the UK working on quantum field theory in curved spacetime: the Cambridge group, 
based around Stephen Hawking and Gary Gibbons; the King's College group based around 
Paul Davies; the Imperial College group, based around Chris Isham and Mike Duff; and 
the Manchester group, consisting of Stuart Dowker and his students. The SLAC preprint 
list would be scanned avidly for anything coming from these groups, but as I was especially 
interested in quantum field theory in topologically non-trivial spacetimes and vacuum en- 
ergy calculations using ^-function methods, the work of Stuart Dowker was always eagerly 
awaited. An especially treasured preprint, that I still use, was Selected Topics in Topology 
and Quantum Field Theory that was based on lectures given by Stuart at Austin, and which 
were never pubhshed; these notes pre-date all of the now standard books and reviews on 
topology and gauge fields. 

Since coming to the UK, Stuart and I have met on many occasions, often at Ph. D. vivas. 
I am sure that many contributors to this volume will have had the "Dowker experience": 
You mention to Stuart that you are working on some calculation to be met with something 
like "Oh yes. That is just a special case of a theorem by Schnekelgreuber^ pubhshed in 
1878." 

^ So far as I am aware I made this name up, but I would not be at all surprised if Stuart was to say 
Scnekelgreuber really existed but was only active between 1829 and 1837, not in 1878. 
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B. Background 



The problem of computing the vacuum, or Casimir, energies in the case of non-smooth 
background potentials has a long history, and in addition has become the focus of much 
recent interest. Of especial interest to the present paper is the case of 5-function background 
potentials. One of the pioneering, and most important, papers on the calculation of vacuum 
energies in 5-function potentials is [1] whose analysis is based on earlier Green function 
calculations of Later work on vacuum energy and other related quantum field theory 



calculations includes 13 



14| . Some motivation for these investigations can be garnered from 



the study of brane-world models stimulated by the Randall- Sundrum scenario [15[ where 
non-smooth backgrounds in the form of 5-function potentials arise in a natural manner that 
is crucial to the features of these models. There are many studies of vacuum energies in 
brane-world spacetimes. (See 16l-l2l| for some of the early investigations.) 

Almost all of the previous studies of vacuum energies in non-smooth backgrounds are 
concerned with non-interacting (apart from interactions with the background) fields. This 
renders the renormalization considerations somewhat simpler than is the case when inter- 
acting fields are involved, particularly beyond one-loop order. The renormalization of 
theory in a (5-function background was considered by Uj but at one-loop order only where 



additional divergences were dealt with by the addition of extra sur 
action. A different interpretation of this procedure was given by 



ace terms to the original 
(see also the later paper 



|l3 |) that is more in keeping with standard renormalization theory, and is the method that 
we will adopt in this paper, to be detailed in Sec. [Ill below. 

One of the key features that we must deal with when proceeding beyond one-loop order 
is the presence of overlapping divergences, and the necessity of showing that they cancel 
(otherwise there will be non-local divergences present that cannot be removed by local 
counterterms) . We will adopt the background field method and dimensional regularization as 
in 22]. A comprehensive review of the method for smooth backgrounds can be found in 23|. 
The general issue of interacting fields on manifolds with boundaries, smooth or not, beyond 
one-loop order appears to have received little attention. (The case of smooth boundaries 
at one- loop order is dealt with in 22|. Tsoupros 24l426| examines smooth spherical cap 
geometries and Haba [^^j considers 5-function backgrounds. The calculations presented 
below do not overlap directly with these references.) 
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The outline of our paper is the following. In Sec. [TTl we discuss the renormalization of 
an interacting scalar field theory with cubic and quartic self-interactions to two-loop order 
in a (5-function potential. Dimensional regularization is used. The background spacetime is 
fiat and assumed to be 4-dimensional with a possible periodic identification of one spatial 
coordinate. The importance of including correct boundary terms, especially a boundary 
interaction in the field (see ( 12. 3 p below), is stressed. Without the proper boundary terms 
the theory will not be renormalizable. Complications in the evaluation of the complete 
two-loop effective action are highlighted, and it is shown how to renormalize the effective 
potential to two-loop order. Necessary two-loop counterterms are found. In Sec. IIIII we 
address the calculation of the one-loop effective potential for a complex scalar field in a 
(5-function potential with one spatial coordinate periodically identified to form a circle. This 
allows the complex field to change by an arbitrary phase around the circle and we calculate 
the effective potential as a function of this phase. A number of asymptotic limits are 
obtained. Sec. HV] summarizes our results and presents a short discussion. Appendix A 
outlines a method for evaluating the Green function in the case of any number of 5-function 
potentials. Appendix B describes some technical details in the evaluation of certain types 
of integrals that arise in the renormalization beyond one-loop order. 



II. RENORMALIZATION 



We consider an interacting scalar field in four spacetime dimensions with a general self- 
interaction. The theory will be regularized using dimensional regularization as in 22|. The 
bare action will be chosen to be 



where 



•S" — S'bulk + S] 



bulk T Obdry ; 



(2.1) 



5, 



bulk 



+ hB<pB + J^9b4>% + \^_^B<P%^ , 



(2.2) 



and 



bdry 



B + 



B + 0-B> ■ 



y=a 



(2.3) 
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The bulk action 5'buik extends over the complete (D + l) dimensional spacetime. The bound- 
ary term Shdrj only extends over the D-dimensional subspace specified hj y = a that gives 
the location of the ^-function. S'bdry could obviously be written in the equivalent form of an 
integral over the {D + l)-dimensional spacetime with the terms in braces in fl2.3p multiply- 
ing the 5-function 6{y — a). All of the coupling constants that occur in the action are bare 
unrenormalized expressions. We might have expected, on dimensional grounds, that a 0^ 
term appeared in S'bdry since its coefficient would be dimensionless for the case of interest 
D = 3. However, this means that the associated 0^ counterterm could not depend on the 
strength of the 5-function (that has units of mass, or inverse length), and because in the 
absence of a 5-function potential there is no need for S'bdry, "we rule out such a term. We 
expect that Sbdry — )■ as f — )■ 0, where v is the renormalized 5-function strength. (This 
conclusion is substantiated in the explicit calculation presented below.) 



In dimensional regularization we require 28|| all renormalized couplings to have the same 



dimensions for all D as they do when D = 3. We write 

D = 3 + e, (2.4) 



and introduce an arbitrary unit of length £. (It is more customary [28|] to use a unit of mass 
fi, but obviously i = l/fi, and we stick with i.) Bare quantities are expressed in terms of 
renormalized ones by 



(pB = 




(2.5) 


ml = 


+ (5m^ , 


(2.6) 


Ab = 


r^(A + 6A) , 


(2.7) 


hs = 


r'/\h + 6h) , 


(2.8) 


9b = 




(2.9) 


Xb = 


tiX + 6X) , 


(2.10) 


Vb = 


V + 6v , 


(2.11) 


IB = 


r^/2(7 + 57), 


(2.12) 


as = 


r'{a + 6a) . 


(2.13) 



We write the field renormalization constant Z as 

Z=1 + 6Z, (2.14) 
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and all counterterms are expanded in powers of h, the loop counting parameter: 



5C = h5C^^^ + ff5C^^^ + ■■■ . , 



(2.15) 



Here C denotes a generic quantity that occurs in fl2.6H2.l4|) . The counterterms 
SC^^\ 5C^'^\ . . . are chosen to contain pole terms in e that cancel the various poles arising in 
the effective action order by order in the loop expansion. 

The effective action can be evaluated in the loop expansion in a suitable form for a 

r 

discussion of renormalization as described by Jackiw | 29 | . Here we follow the formalism 
used to consider renormalization in curved spacetime 22|, |30[ . The only added feature here 
is the presence of 5'bdry This affects the Green function, as well as the interaction vertices 



beyond that occurring for the bulk theory. The basic step is to write 



(P = ^ + n}i^ij , 



(2.16) 



with Lf the background field. The results (12.5112.14]) are used in (12. 2 p and (12. 3p . all coun- 
terterms are expressed as in (I2.15P and the result expanded out to a consistent order in h. 
Standard functional methods then generate the loop expansion of the effective action. The 
result can be expressed as 



m = S[^] + ^trln(£2A) + h/l- exp (-^^int 



(2.17) 



where A is the operator that defines the Green function following from the second functional 
derivative of the action functional with respect to the field evaluated at the background 
field. (See (12.181) and (I2.19P below.) The terms of cubic and higher orders in the field, as 
well as the counterterms, are defined as the interaction part of the action S'int and treated 
perturbatively. The angular brackets in (I2.17P denote an evaluation using Wick's theorem 
with only one-particle irreducible graphs kept. We define 



{(t){x)(f){x')) = Gy{x,x') , 



(2.18) 



where 



-□^ + 171^ + g(p + + v6{y - a) 



Gv{x, x') = 6{x, x') 



We can express S'int in terms of h by defining 



h 



(2.19) 



(2.20) 



In our case, 

Ai = j d^+'x ^ {g + A^(x)} , (2.21) 

J d'^x • (2.22) 

y=a 

The only difference with what happens in the absence of a (5-function potential is the presence 
of the last term in ^42, and an alteration of the Green function. (Of course the alteration of 
the Green function affects the evaluation of the effective action considerably.) 

Expanding the effective action to two-loop order results in (see 3l| for a review) 

r[^] = S[^] + hrWK) + ({A,) - hAl)] + ■ ■ • 



where 



S[^] + ^tr In (£2 A) + h^T^^^ + ■ ■ ■ , (2.23) 



r(2) = [ d''+'x\^Gl{x,x) + l-6m"^'^G,{x,x) 



^6g'^'^^{x)G,{x, x) + ^6X^'^f{x)G,{x, x)} 



+ 2 

+ j d^x ]^5v^^^G^{x,x) 

y=a 



- ^r^'J d''+'xd''+'x'[g + X^ix)][g + X^ix')]Glix,x') , (2.24) 

gives the two-loop contribution to the effective action. The next two subsections will examine 
the one- and two-loop divergences of the effective action that we have obtained. 



A. One- loop effective action 

The one-loop part of the effective action fl2.23p involves 

r« = -trln(£2A) . (2.25) 

In order to deal with this expression using dimensional regularization, we can differentiate 
with respect to and use the fact that the inverse of A is the Green function as defined 
by <^M- We find 

^.r(i) = i / d^+ia;G.(x,x). (2.26) 



It is sufficient to compute G^i^x, x) with constant background fields, since using the derivative 
expansion and power counting shows that there can be no pole terms that involve derivatives 
of the background field at one-loop order. (This reflects the well-known result that field 
renormalization is ffist required at two-loop order. See IsL] for a proof of this using the 
derivative expansion of the effective action.) We can then use the result in (1A3P with m? 
replaced by 

^ ~ (2.27) 



ip is now regarded as constant. Use of (1A3I) gives us 

L/2 



dm? 2 ^ j ^ 

-L/2 



(2^ 



Gv{p;y,y) 



(2.28) 



where Vd = J d^x is the volume associated with x_|_. If we use the result in (]A16P the 
integration over y may be performed with the result 



d 



dm' 



d^p 

(2^ 



L 

2uJr 



coth(La;p/2) 



V [Lujp + sinh(Lco'p)] 



8a;3sinh^(Lwp/2) 



(2.29) 



Note that because of the replacement of m^ with in ( I2.27p . we should use Up = 
(p2 + M^)^/^ here. 

If we are only interested in the pole part of T^^^ to discuss the renormalization, then we 
may study the large p behaviour of the integrand. If we denote the pole part of T^^^ by 
PP(r(i)), then it follows that 



d 



dm^ 



pp (r<") 



d^p 

(2^ 

V 



2w„ 



1 



2Ur. 



(2.30) 



Terms that have been dropped here decay exponentially for large p. Although they con- 
tribute to the finite part of T^^\ they make no contribution to the poles. Using the definitions 
in f IBip and f IBSp . and using the recursion relation f lB4p shows that 



d 



dm'- 



.pp (pp)) 



ivo PP (^/(i) - j/(3) 

+y/(4)-^/(5) + .. 



(2.31) 



Because of the simple dependence of on m^, using ( ]B2p results in 



pp (r«) = -Av-^(4^)-^/2 pp 



L r(-l/2-D/2) 

2 r(i/2) ^ ^ 



t;r(l/2-D/2),,^2,D/2-i/2 , ^'r(l-Z}/2),,^2,D/2-i 



+ 



4 r(3/2) ' ' 8 r(2) 

T;3r(3/2-D/2),,^2,^/2_3/2 



+ 



(2.32) 



16 r(5/2) 

The next term to that indicated will be finite as D — ?■ 3 and therefore contains no pole. By 
taking = 3 + e and expanding about e = we find 



PP (r«) 



^3 

967r2e 



(3LM^ + ^vM'^ - 



(2.33) 



Given the form of the counterterms in f l2.5H2.15|l . with bZ'^-^^ = and tp constant, we have 
the one-loop counterterm part of the action as 



(2.34) 



Note that J d^x = V^L identifies the bulk part of the action and the L-independent part 
identifies the boundary part. Using the expression for in ( I2.27P in the pole part of the 
one- loop effective action in f l2.33p . it is easy to show that the resulting effective action is 
finite as e — ?■ if we choose the counterterms to be 

,4 



m 



gm? 
1 

3A^ 
\v 

V 

9 



(Am^ + g') 



[m 



(2.35) 

(2.36) 

(2.37) 

(2.38) 

(2.39) 

(2.40) 
(2.41) 

(2.42) 



The counterterms that enter the boundary part of the action all vanish as — )■ as would 
be expected. The counterterms that enter the bulk part of the action agree with the more 
general case considered in Ref. I22I. 

B. Two- loop effective action 

The complete two- loop effective action was given in (12 .24^ . It has not been possible to 
evaluate even the pole part of this expression for the general scalar field theory that we have 
been considering. The difficulty with performing a complete renormalization calculation 
concerns the extraction of the complete pole part of Gl{x,x'), a problem that we were not 
able to overcome due to the complicated nature of the Green function. Instead we will 
examine the simpler, but still significant, task of renormalizing the vacuum energy for a 
X(f)^ theory obtained by setting h = g = and taking the background scalar field to vanish 

= 0). With this simplification, we find 

j d^x^6v^^^G^{x,x) . (2.43) 

y=a 

We have already described how to extract the pole part of G^{x,x) in Sec. Ill A[ At two- 
loop order we have the one-loop counterterms Sm?^^^ and 5v^^^ multiplying Gv{x, x), resulting 
in poles that involve the finite part of Gv{x,x). Such poles, should they not cancel, will 
involve complicated non-local expressions and would render the theory non-renormalizable 
with local counterterms. A necessary part of the analysis will be to show that all such 
non-local poles cancel between the three terms of F*-^^ given in (12.431) . 

We have given Gv{x,x) in (]A17|) with Gv{p;y,y) given in (1A16P for the case of peri- 



odic boundary conditions. Because our main focus is on how the presence of a 5-function 
potential changes the renormalization calculation, in order not to complicate the analysis 
unnecessarily, we will let L — t- 00. (The presence of a finite L leads to more complicated 



non-local divergences as known in the absence of 5-function potentials from 32|, |33| . It is 
possible to show that in this more complicated setting all non-local divergences still cancel 
as we describe below for infinite L, although we omit the details of this here for simplicity.) 
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With the Umit L — )■ oo taken in ( ]A16P we have 



This is the resuh that we would have found had we simply chosen not to adopt periodic 
boundary conditions over a finite interval and taken the whole real line instead. 

The evaluation of the Green function expressions that enter the two-loop part of the 
effective action in fl2.43p is described in Appendix |Bl The first two terms, that contain an 
integral over y as well as over x_l, give rise to bulk counterterms (ones that multiply LV^) as 
well as surface counterterms (ones that multiply V3 but that are independent of L). Because 
the one-loop counterterms Sm?^^^ and Sv^^^ multiply Gv{x,x), we must include the finite 
part of G^{x,x) in order to calculate the complete pole terms of the effective action. From 
(IB2D . with L> = 3 + e, we have 



2 2 
m m 



1(1) = ^ + 



An^e 87r2 



+ ■ ■ ■ , (2.45) 



/(3) 



+ • • • , (2.47) 



27r2g 47r2 

upon expansion about e = 0. (The ellipsis indicates terms of order e that cannot lead to 
poles in the effective action, although they may contribute to the finite part.) 

If we temporarily ignore the surface term in T^'^\ and use the expression found earlier for 
Sm^^^^ given in 02.371) . it is possible to show that 



Xv^ / Xvm'^ Xv^ 
1024^ ^ V 512^ ~ 10247r4e / V An 

XV^ /• .n 1 



KJv) + r' [ d^x-Sv^^^GJx,x) (2.48) 
512n^e ^ ' J 2 ^ ' ' ^ ' 

y=a 

contains the pole part of the two-loop effective action. The bulk term, the first term propor- 
tional to L, involves only a local divergence that can be dealt with by a local counterterm. In 
this bulk term, all of the dependence on 7 and ln{m'^P) that occur at intermediate stages of 
the calculation has cancelled, but there are observed to be divergences that involve ln{m'^i'^) 
in the surface term. In addition, there is a non-local divergence proportional to the integral 
K4{v) present. 
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The surface contribution (the last term on the right hand side of fl2.48p ) involves 
1 





Xvm? 


5127r4e2 


5127r3e 


Xv^ 


/ Xv^ 


' 5127r4e2 + 


Vl0247r4e 


Xv^ 




+ 512vr2e^^' 


» + ■•• , 





Xvm? \ ^ ( vr?"'^ 

1^ ■ A 

4-71 



(2.49) 

if we use fl2.40p for and flB6l) . If we use the result found in fl2.49p back in (12.481) it can 
be seen that the non-local divergences that involve K4^{v) cancel, and that the dependence 
on ln(m2£2) also cancels. We are left with the simple expression 

/r,N , Am^ Xv^ Xvm? 

as the pole part of r*^^) ^he first term on the right hand side represents the bulk divergence, 
present even when there is no 5-function potential, and the other two terms represent the 
poles on the surface. Note that the presence of the contribution from the one-loop sur- 
face counterterm, that involved the renormalized strength of the 5-function, was crucial in 
obtaining a cancellation of the non-local divergences. 

The divergences present in (I2.50p can be seen to be removed by the choice of two-loop 
vacuum energy counterterm 



and the renormalized surface density 



2567r4e2 V 4 

This presents a complete proof that to two-loop order the vacuum energy density is renor- 
malizable, and we have computed the necessary counterterms to do this. 

III. VACUUM ENERGY 

In this section we will compute the vacuum energy density to one-loop order for a complex 
scalar field $ in a 5-function potential. We will allow the boundary condition on the field 
to depend on an arbitrary change of phase around the circle: 

$(-L/2) = e^^™$(L/2) . (3.1) 
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Here < a < 1 is the arbitrary phase factor. The case a = corresponds to periodic 
boundary conditions, while a = 1/2 corresponds to antiperiodic boundary conditions. We 
will look at how a enters the expression for the vacuum energy density. 

Similar calculations without a 5-function potential present were performed some time 
ago. Of special interest is the paper of Ford [3j] , where it was remarked that the situation 
is the same as coupling a scalar field to a constant gauge potential, with the case of a = 



corresponding to a vanishing gauge field unstable 
to non-Abelian gauge theories 35|], and Hosotani 



br spinor fields. This was later generalized 



36j showed that by computing an effective 



potential as a function of a it was possible to break non-Abelian gauge symmetries. More 



recent work [SH] has studied this mechanism in the brane-world setting where 5-function 
potentials arise naturally. 

The generalization of (12.261) to the complex scalar field case is obtained simply by multi- 
plying by two (to account for the fact that a complex field can be represented by two real 
fields) and changing the boundary conditions on the Green function to reflect (13. ID : 

L/2 



d 



dm- 



tw = V, 



D 



dy 



-L/2 



(2^ 



Gv{p;y,y) 



From (1A12P we find 



Gv{p;y,y) = Go{p]y,y) 
where Go{p; y, y') is given by (]A13|) with 



2n 



^3 



vGo{p;y,a)Go{p; a,y) 
1 + vG(,{p; a, a) 



(3.2) 



(3.3) 



(3.4) 



and Up=p'^ + rn? as before. From (1A13P it is easy to see that 



L/2 

f °° 1 

/ dyGo{p;y,a)Go{p;a,y) = ^ -(a| + Wp)" 



-L/2 



d 



dm? 



oo _ 



^2 I , ,2\~1 



(3.5) 



This leads to 



r(i) 



surf ' 



(3.6) 
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where 



(3.7) 



j=-oo 



is the term present in the case of f = 0, corresponding to the absence of a (5-function 
potential, and 



^ surf - 



(2vr) 



D 



In 



2)-l 



j=-oo 



(3i 



is the 'surface' contribution to the effective action coming from the 5-function potentiaL 
The sum over j in fl3.8p is easily evaluated by contour integral methods [38| to give 



2\-l 



sinh(Lci;p 



j=-oo 



2ujp cosh.{Lup) — cos(27rQ;) 



(3.9) 



The pole part of T^Ji^^^ as D — )■ 3 can be analyzed as described in Sec. [TTl above. 

The r^^Q term can be related to the function evaluated first by Ford [34]. (See [311 ^ 
textbook treatment.) If we define 



F{\-a,b)= J2 [{n + af + b'']' 



(3.10) 



then after suitable analytic continuation 



where 



fx{a, b) = 4sin(7rA) j dx {x^ - fo')"^?? [e2'^(-+-) _ i] ' 
b 

It is straightforward to show that 



(3.11) 



(3.12) 



v=0 



.^..-(4.)-r(-f)(|)%(-f:.,- 



27r 



(3.13) 



and that poles come only from the first term on the right hand side of (13. lip when used in 

(EH. 

A full renormalization calculation can be performed as we have described in Sec. HTl 
However, our aim in this section is to study the role that a and v play in the expression 
for the vacuum energy density. To do this we will simply compare the vacuum energy for 
general a with that for a = 0. (A detailed analysis of the a = case was given in [l^ 
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and will not be repeated here.) Because T^^^ = V^Lp with p the vacuum energy density, or 
effective potential, we have 

p = p(a)-p(0) (3.14) 
as the difference between the a 7^ and a = cases. A bit of calculation shows that 



00 



mL 

/ y + vL/2 - 2ye-y cos(27ra) + - vL/2)e-^y \ 

1 (1 - e-y) [y + vL/2 - {y - vL/2)e-y] / ' ^ ' > 

This same result can be found by the more tedious and lengthy process of expanding about 

the pole at D = 3, removing the poles with counterterms as described in Sec. HIl or else by 



13|. 



adopting (^-function regularization where no poles occur 

If we view p as a function of a, it is easy to show from fl3.15p that for < a < 1, p has a 
global minimum at a = (or a = 1) and a global maximum at a = 1/2. Thus, the case of 
antiperiodic boundary conditions leads to the maximum energy density, a conclusion that is 



the same as in the absence of a 5-function potential [SJ, l39|| although the actual expressions 
for the energy density differ of course. In Fig. [1] we show the result of evaluating f l3.15p 
numerically for different values of mL, but with vL kept fixed. The vacuum energ y d ecays 



35|. 



exponentially with mL exactly as in the case with no 5-function potential present 

We can also study what happens if the strength of the 5-function potential v is varied. 
For simplicity we will set m = and a = 1/2. The result is plotted in Fig. [2] as the solid 
curve. If we take fL ^ 1, it is possible to obtain the following asymptotic expansion for p 
in f l3.15p for m = but a general: 

vr^ sin^fvra) r 4 



_ vr sm (TTa) f 4 ^ / ^ 



2 

[86 + 3tt^ + 2(8 - TT^) cos(27ra) + (10 - vr^) cos(47ra)] 



3v^L^ 

+ ■■■}. (3.16) 

For the case of a = 1/2 this is plotted as the dotted line in Fig. [2j 

We can also analyze the case where vL is small and obtain a reliable asymptotic expansion. 
The analysis is reasonably involved, so we will simply quote the result. If vL/2 ^ 1 we find 
from fl3.15p that (again taking m = and a = 1/2 as an example) 

P^^-^vL + livLf^ - ^{vLf + . . . , (3.17) 
15 




FIG. 1: This shows plots of p as a function of a, the parameter that determines the boundary 
conditions, for three different values of mL. The top curve is the result for m = 0. The middle 
and lower curves are the results for mL = 1 and mL = 3 respectively. In all cases the maximum 
occurs at a = 1/2 as the analytic proof described in the text shows. The energy density decays 
exponentially with mL. 

where (3 ~ 0.0288734 is the result of a numerical evaluation of a simple integral of the expo- 
nential integral form. Contributions of the next order in fl3.17p can also be found and involve 
non-analytical terms of order (fL)^/^ and (fL)^ln(fL). The fact that p is not analytical at 
f = is what complicates attempts to calculate the asymptotic expansion. Similar non- 



analyticity has been seen in the earlier calculations of [lOj in similar settings and indicates 
the futility of trying to use a naive perturbative approach about v = Q (corresponding to 
treating the 5-function potential as a perturbative interaction.) To demonstrate the utility of 
fl3.17p we plot the approximation (shown as the dotted line) against a numerical evaluation 
of the exact result (shown as the solid line) in Fig. [3l 
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5 10 15 20 25 30 

vL /2 

FIG. 2: The solid line shows the result of 27r^L^p as a function of vL/2 found from (I3.15P in the 
case m = and a = 1/2. The dotted line shows the same result found using the asymptotic 
expansion in ()3.16p . The agreement between the analytic and numerical results become very good 
once vL is sufficiently large. 

IV. DISCUSSION AND CONCLUSIONS 

We have considered the case of an interacting field theory in a non-simply connected 
spacetime in the presence of a (5-function potential. The necessity for a proper inclusion of 
boundary interaction terms for renormalizability was discussed. The required counterterms 
were computed, using dimensional regularization, to one-loop order, and some consideration 
was given to the two-loop counterterms. It was shown why the complete two-loop calculation 
was difficult. We then showed how to obtain the effective potential for a complex scalar 
field with general boundary conditions around a compact spatial dimension, generalizing 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 

vL /2 

FIG. 3: The solid line shows the result of 27r^L^^as a function of vL/2 found from (I3.15P in the case 
m = and a = 1/2. The dotted line shows the same result found using the asymptotic expansion 
in ()3.17p . Note that the range is shown over 6 — 8.5 rather than extending to to exaggerate the 
difference between the two curves. As vL decreases, the agreement becomes excellent. 

earlier such studies. A number of approximations were obtained using analytical methods 
and compared with a numerical evaluation of the exact result for the effective potential in 
certain cases. In particular, the case of a weakly coupled 5-function potential was shown to 
result in a non-analytic expression for the vacuum energy that will not show up if a normal 
weak-field perturbative approach is used. 

There are a number of future directions that are worthy of attention. The first is to find 
a method to obtain the complete set of boundary counterterms to two-loop order, and if 
possible to proceed beyond two-loop order. In particular it would be of interest to see if 
the non-analyticity seen in the 5-function vacuum energy could affect the renormalization 
procedure at higher orders. It would also be of interest to examine how the analysis pre- 
sented here is modified if there is more than one 5-function present, or if the 5-functions 
are more than one-dimensional. The case of spherical 5-functions relevant for spherical or 
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cylindrical boundary problems is also of some interest. Having a more complete analysis of 
the counterterms would enable a renormalization group study to be performed and further 
illustrate the role of the surface divergences. Finally we mention that it is of interest to 
examine the complete stress-energy tensor for the interacting case, complementing the free 
field cases that have been done. Some work has been done on this 43 1. 
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Appendix A: Green functions 

In this appendix we describe a method for calculating the Feynman Green function in the 
presence of one- dimensional 5-function potentials that generalizes ^|. An earlier evaluation 



of the one-dimensional Green function us,ng an entirely different approach was given by |40|, 

and by [2| for scalars and spinors. Take a (D + l)-dimensional spacetime with spacetime 
coordinates = (x^,?/) and adopt a Euclidean metric. Here y is used to distinguish the 
coordinate that enters the 5-functions. Take the potential to be 

n 

V{y) = J2v,S{y,ai) , (Al) 
1=1 

where the Vi are constants and give the locations of the (5-function singularities. In the 
case where y runs over the finite range [—L/2, L/2] with y = —L/2 identified with y = L/2, 
we assume that —L/2 < < L/2 for alH = 1, . . . , n. 

We now wish to solve for the Green function Gy{x, x') defined as the fundamental solution 

to 

(-□ + m^ + V{y)) G^{x, x') = S{x, x') . (A2) 

Write 

G.{x,x') = I ^ e^^<--'-'^^G.{p;y,y') . (A3) 
It then follows from ( IA2|) that 



Oy2 +^P + viy) 



G,{p;y,y')=6{y,y') , (A4) 
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where 



•p = {p^ + m^) 



(A5) 



With V{y) given as in flAip we can use the 5-functions to set y = ai in G^ip; y, y'), thereby 
obtaining 

92 



dy^ + 



Gv{p;y,y') = S{y,y') -'^Vi6{y,ai)Gy{p; 



(A6) 



i=l 



Now define Go{p;y,y') to be the solution to (1A6P with t>j = 0; that is, for no (5-functions 
present. This means that 



Go{p;y,y') = S{y,y') . 



(A7) 



We can set y' = ai in (IA7p and then use this to ehminate 6{y, ai) in ( \A6\i . Rearranging the 
resuh gives 



dy^ + 



Gv{p;y,y') + Go{p;y,ai)G^{p; ai,y' 

= S{y,y'). 



i=l 



(A8) 



Since the solution to (]A7|) should be unique (given the boundary conditions) we may identify 
the expression in braces in flASP with Go{p; y,y'). This gives us 



Gvip;y,y') + ^ViGo{p;y,ai)G^{p; ai,y') = Go{p;y,y') 



(A9) 



1=1 



If we set y = a,- in ( 1A9I) we find 



^ {5ij + ViGoip; ai)} Gy{p; ai, y') = Gq{p] aj, y') 



(AlO) 



i=l 



This gives us a set of equations that determines G^{p; ai,y') that occurs in (\A9\i in terms 
of Go{p;y,y'). Thus we obtain Gy{p;y,y') from a knowledge of the Green function in the 
absence of a (5-function potential. 

In the simplest case of a single 5-function potential (n = 1 in (lAlOl) ) it is easy to see that 

Go{p;a,y') 



This gives us 



Gv{p;a,y') 



Gyip;y,y') = Goip;y,y 



1 + vGo{p; a, a) 

V Goip;y,a)Go{p; a,y') 
1 + vGo{p; a, a) 



(All) 



(A12) 
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if we use ( ]A9p . The case of more than one (5-function potential can be dealt with in a 
similar manner, although of course the details become more involved. For exarnple, it is 
straightforward to recover the Green function for two 6 functions used by Milton 10|. 

For the present paper we are concerned with a single 5-function, and assume —L/2 < y < 
L/2 with the endpoints identified. For the case of periodic boundary conditions we have 



oo ^ 



where 



2nj 



(A13) 



(A14) 



Note that Go{p; y, y') can only depend on \y — y'\, a result that is easily seen by relabelling 
i to — ? in the sum. The sum over j can be computed using contour integral methods to 



give 



4l| 



Goip;y,y') 



1 cosh {^Up - ujp\y - y'\) 



(A15) 



2up sinh (f^p) 

This expression is sufficient to determine Gy{p;y,y') and hence the full Green function in 
the presence of a 5-function potential. Note that as L — )■ oo we recover the flat spacetime 
result of Go{p] y, y') = e-^^\y~y'\ / {2up). 

We are mainly concerned with the coincidence limit of the Green function. From (1A12P 
and ( 1A15I) we find 

V Gl{p;y,a) 



Gv{p;y,y) = Go{p;y,y) 



1 , fL 

coth — Wr, I — 

2uj„ V 2 ^ 



1 + vGq{p] a, a) 
V cosh^ [("I - \y - a\) Up] 



H smh'i^up) [l + ^coth(fc.,) 
The coincidence limit of the Green function is then obtained from ( ]A3I) . 



(A16) 



with ( ]A16|) used on the right hand side. 



(27r) 



D 



Gv{p;y,y) 



(A17) 



Appendix B: Evaluation of some integrals 



We define 



1(a) 



d^p 
(2^ 



(Bl) 
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where ujp = ijP' + w?)^/"^. Making use of the standard integral representation for the F- 
function 4l|] it is easy to show that 



/(a) = (47r) 



-D/2 



r(f) 



m 



2^{D-a)/2 



(B2) 



This is a standard result of dimensional regularization 42 1. 
We also define 



KJv) 



{271 



1 + 



2Ur. 



It is simple to show that Kn{v) satisfies the recursion relation 



KJv) 



I{n) - -K„ 



(B3) 



(B4) 



where I{n) is defined as in fIBip . We are concerned with the case — )■ 3 in this paper, and 
simple power counting shows that Kn{v) is finite as D — ^ 3 for n > 4. The recursion relation 
( ]B4[) enables us to isolate the divergent parts of Kn{v) in terms of the simpler integral given 
in (EI]) and flB2l) . 

With these preliminaries over, we can now evaluate the divergent parts of the Green 
function expressions that enter the two-loop vacuum energy. First of all, in the large L 
limit, from (1A16|) and (1A17P we have 

L/2 



dy Gj;{x, x) 



-L/2 



d^p 

(2^ 



L 



4u;2 



1 + 



2UJr 



dye 



-2\y—a\ujp 



d^p 

j5 



4 J i2n 



-1 



2 ^ ' 



(B5) 



In the last line we have used (1B4I) . The result in (IB 5 1) is exact. The first two terms contain 
poles as — ^ 3 and the last term is finite. 
We also need 



y=a 



d^p 

(2^ 



1 V „, 
2^ ~ i'^*' 



1 + 



2UJr 



-1 



1/(1) - "^K.iv) 

- \i{2) ^ 



(B6) 
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Any poles can come only from the first three terms on the right hand side as K^lv) is finite 
as D ^ 3. . 

L/2 

Finally, we need J dy Gl{x, x). We can use (IA16P and flA17P to show in the limit of 



L — )■ cxD that 



L/2 



L/2 



dyGl{x,x) 



-L/2 



dy 



d'^p d^q ( 1 



-L/2 



(27r)^ (27r)^ UwpW, 



1 + 



2UJr. 



-2\y-a\ujq 



2 / \ -1 

2 2 / . ^ 



1 + 



2Wr 



-2|j/-a|(t<^p+(^9) 



(B7) 



with 



J = ^ 



16 J {2ti)D {2ti)'^ 'p 



2uj, 



2uj.„ 



(B8) 



The only complication is the evaluation of J, because the double integral does not factorize. 
If we consider the integral over q first it is easily seen that 



d^q 



{2n 



2Ur. 



V 



2 ( 

v\-^ f d^q ^1 



[^P 2 



{2ny 



(B9) 



If we use flB9P in fIBSP , we find that 



J = ll{l)Ks{v) - j-K,iv) [K,{v) + K,{~v) - 2/(1)] + J, 



(BIO) 



where 



with 



2\ -1 



Flip) 



d^q 
(2^ 



(Bll) 



(B12) 
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We only require the pole part of J for the divergent part of the effective action. It is easy 
to show that 



with 



F2{p) 



(27r)^ " 

F2{p) is finite as D — )■ 3, and in setting D = 3 we find 

1 



(B13) 
(B14) 



m _. 



m 

+ —UJ, 



7 CO. 



-1 



27r2 P 



27r2 P 

Using flBT3|) and flBTSl) in flBTTjl results in 



1/2 



In 



+ 



1/2- 



(B15) 



2 3 

647r ^ ' 327r2 ^ ^ 
[K,{v) + K,{-v)] 



128 

'2567r 



(27r 



X In 



+ 



m 





2567r2 










(4- 


s 1/2 


■0 



1 



1/2 



(B16) 



This result is exact. We are only after the pole part of J, and hence Ji, so we can drop 
all terms that are finite as — t- 3. Of the expressions /(l),/(2) and /(3), only /(I) and 
/(3) contain poles as — > 3. K/^{-^v) is finite as D — i- 3. The only remaining question 
concerns the last term in (IB16P as -D — ?■ 3. To see the pole structure of this integral, expand 
the integrand in powers of ujp. Terms that fall off with Up faster than uj~'^ will converge as 
D -> 3. We find 



m 



{2n 



V 



d^p 



1/2 



In 



/ , ,2 \ 1/2- 

m \ m'^ 



(B17) 



(27r)^ ^ P 

where terms that are finite as D — > 3 have been dropped. From (IBip it is easy to see that 

d^p 



(B18) 
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Because /(2) is finite (after regularization) as D — )■ 3, so is /'(2). Tliis means tliat (1B17P 
is finite after regularization as -D — t- 3. We conclude tliat tlie last term of (IB16P does not 
contain any pole terms. 

Retaining only terms that can contain poles as D — )■ 3, we find from flBlOp that 

2 2 4 

?? 7? Tnil 11 

J = -mm - -m - ^m) + ^mK,{v) + ■■■ . (bi9) 

Liberal use of the recursion relation (]B4p has been made here. It is noteworthy that all of 
the terms that involve K{—v) at intermediate stages of the calculation have cancelled. The 
result for f lB7D becomes 

dyGl{x,x) = ^/^(l) - ^/(l)/(3) + ^1(1)7^4^ + ^/(2)/(3) 

2 4 
1) Tfil) 1) 

Again terms that are finite as — t- 3 have been dropped. The presence of non-local pole 
terms coming from 1(1) and 1(3) multiplying Ki{v) can be noted. Such terms must cancel 
for renormalizability with local counterterms, and it is shown in Sec. Ill Bl that this is the 
case. 



L/2 



-L/2 
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